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q . Abstract: A novel and simple approach to optical wavelength measurement is 

jy! ' presented in this paper. The working principle is demonstrated using a tunable 

Q . waveguide micro ring resonator and single photodiode. The initial calibration is 

55 ' done with a set of known wavelengths and resonator tunings. The combined spectral 

£>->, sensitivity function of the resonator and photodiode at each tuning voltage was 

rG ' modeled by a neural network. For determining the unknown wavelengths, the 

resonator was tuned with a set of heating voltages and the corresponding photodiode 
signals are collected. The unknown wavelength was estimated, based on the collected 
photodiode signals, the calibrated neural networks, and an optimization algorithm. 
The wavelength estimate method provides a high spectral precision of about 8 pm 
(5T0 -6 at 1550 nm) in the wavelength range between 1549 nm to 1553 nm. A 
^) ' higher precision of 5 pm (3 TO -6 ) is achieved in the range between 1550.3 nm to 

\Q , 1550.8 nm, which is a factor of five improved compared to a simple lookup of data. 

The importance of our approach is that it strongly simplifies the optical system and 
f^) , enables optical integration. The approach is also of general importance, because it 

CO ' may be applicable to all wavelength monitoring devices which show an adjustable 

wavelength response. 
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1. Introduction 

Measuring and controlling the wavelength of lasers is essential to a vast number of applications. Examples 
range from multichannel wavelength division-multiplexing (WDM) JT|, optical communication |2), linear 
and nonlinear spectroscopic applications ||3}|6] , to laser b ased metrology f 7) . It is desirable to combine high 
precision and simplicity with small size and the option of integration. 

State-of-the-art wavemeters, such as double-folded Michelson interferometers JS], can readily provide 
a high spectral resolution of better than 10~ 6 . However, this comes at the cost of a fairly large instrument 
size (in the order of 10 6 wavelengths, i.e., typically a meter). 

There are also simplified wavemeters which are suitable for miniaturization. Their working principle 
is based on two (or several) channels in which the transmission vs. wavelength is different (5). If the 
mapping between the transmission of the channels and the wavelength is bijective (i.e., the transmission 
curve of each channel is strictly increasing or decreasing) , the normalized transmission ratio of the channels 
can be calibrated in the form of a look-up table (LUT). Miniaturization can be achieved by using the 
transmission functions of integrated optical wavelength filters IfTOl , such as thin-film interference filters, 
photonic crystal waveguides ifTTl . or multimode interference couplers |[T2l . This also allows the laser and 
its wavelength monitor to be integrated into a single device |fl3ll . However, such wavemeters deliver only 
moderate precision and the operational range is confined to the region where the spectral sensitivities of all 
channels are either strictly increasing or decreasing. 

In this paper, we present a new method for wavelength and power estimation, based on the calibrated 
transmission spectra of a micro-ring resonator (MRR), providing high precision and extended measure- 
ment range suitable for integrated optics. A wavemeter was constructed in an integrated optical design for 
verification of the method, which simply consists of a single, on-chip tunable MRR as an optical transmis- 
sion filter and a single photodiode. The MRR used here is fabricated from SisN/t/SiCh with a box-shaped 
waveguide cross section (TriPleX™) fl4l . The wavelength and power of the incident light is estimated by 
analyzing the photodiode signals with neural networks and a nonlinear optimization algorithm. 

The contribution of our work lies in the method for high-precision, range-extended wavelength and 
power estimation to simplify the optical requirements. Furthermore, with the proposed wavelength esti- 
mate algorithm, the spectral sensitivity of the MRR is not necessarily required to be strictly increasing or 



decreasing in the operational range. Therefore, the operational range is extended and more tunings can 
be applied for wavelength estimate (i.e., more heating voltages can be applied to the heater of the MRR), 
which results in a higher precision in the measurement. An advantage of this wavelength estimation method 
is that monotonicity is not required, allowing one to consider a larger number of physical implementations, 
including devices with randomly varying transmission. Thus, one may consider the transmission spectrum 
of a powder lfT31 . or a multi-mode fiber lIT6l . for example, as long as the transmission function is repro- 
ducible. 

2. Principle of wavelength and power estimation 

Fig. Q] shows the schematic of the wavemeter. Laser light with an unknown wavelength X x and unknown 
power P x passes through a tunable optical filter. The transmission of the filter vs. wavelength, /(A , v k ), can 
be modified with an external control parameter, vj. In our case, the filter is a single MRR equipped with an 
electric heater (see Fig J2}. The control parameter is the voltage applied to the heater, which tunes the optical 
resonator length of the MRR and thereby modifies its transmission function. To estimate the wavelength 
of the input light field, a number of Af different voltages are applied, where we define v k as the k th voltage. 
The transmitted light illuminates a photodiode, which yields a set of measurement signals y^, in our case 
voltages. The sensitivity of the photodiode is denoted as d(X), which is a function of the wavelength A. The 
measurement signal y k is given by the product of the incident laser power, the transmission of the MRR, 
and the sensitivity of the photodiode. When including measurement noise, r\ k , the measurement signal can 
be expressed as 

Jk = P x f(X x ,v k )d(X x ) + ri k = P X S(X X , v k ) + ri k . (1) 

Here we have abbreviated the product of the MRR transmission and photodiode sensitivity as the combined 
spectral sensitivity function S(X X , v k ) = f(X x , v k )d(X x ). 
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Fig. 1. General scheme of the wavemeter. Laser light with wavelength X x and power P x 
passes the tunable color filter, i.e., a micro-ring resonator (MRR) in this paper. The spectral 
transmission function of the MRR can be changed by applying a heater voltage v k . The 
transmitted light is collected by a photodiode, which yields a measurement signal v^. 

To obtain, from measurements of y k , the unknown wavelength, X x , and power, P x , of incident light, we 
proceed in two steps: 

1. Device calibration 

The first step is to obtain information on the spectral sensitivity function. In principle, S(X,v) may 
be derived via physical modeling of the MRR transmission and the photodiode response. However, 
in practice, the highest accuracy is achieved via calibration. Calibration requires injecting light at a 
known power at a number of known input wavelengths, A, into the MRR. For each wavelength, A, 
the heating voltage, v k , of the MRR is tuned to N values, i.e., k = 1, — ,iV with N the number of 
tunings; and the corresponding photodiode signals, y k , (k~ I,- ■ ■ ,N) are recorded for each pair of A 
and Vk. 

Fig- El shows typical sets of photodiode measurement data. Each data set is recorded at a different 
heater voltage v k . An approximation for S(X,v) can then be obtained by fitting an analytical function 
to the data ifTTl . such as a polynomial, a spline or a neural network (NN). A 2-layer neural network 
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Fig. 2. Microscope picture of the waveguide chip (left) with micro-ring resonators (MRRs) 
and a schematic drawing (right) of the same chip. "R" denotes the MRR with the heater 
on top (gray). "E" denotes electrical contacts (gold colored) and the black lines represent 
waveguides. The fiber from the tunable laser was connected to the MRR from "IN". "PD" 
denotes the position of the photodiode. 
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Fig. 3. Spectral sensitivity function S(A , v) at different v. "+": calibration data; lines: curves 
fitted by neural networks. For clarity, curves are shown only for five heating voltages 
v =0 V, 2.5 V, 5 V, 7.5 V and 10 V. 



is able to model a broad range of nonlinearities 1184-1201 and, from a practical point of view, can be 
implemented and trained very conveniently, e.g., with a neural network toolbox lETI . Therefore a 
2-layer neural network with Q neurons in the first layer and one in the second layer was chosen for 
our work. At a certain control voltage v k , the output y k (X) of the neural network is given as 



y k (X)=P-S(X,v k ) =PY, w iiktanh(w 2 ikX+s uk ) + s 2 ik- 



(2) 



;=1 



Eq. (O models the spectral sensitivity curve as a superposition of step response functions (neurons), 
in this case chosen as tangent hyperbolic functions, wuk an d w 2ik are the input and output weights of 
the neural network, respectively; sijj and 52,1 are biases on the input and output neurons, respectively. 
To obtain a good compromise between a best fit and a minimum number of fit parameters, the number 
of neurons, Q, should be selected according to the required fitting accuracy of the neural network. 

Wavelength and power estimation 

In the second step, light with an unknown wavelength, X x , and power, P x , is sent into to the waveme- 
ter. A set of N photodiode signals, y k , is recorded via applying N different heater voltages, v k 
(k= 1, • • • ,N). The heater voltages can be equally spaced (as in the experiments described in this 



paper), randomly selected or designed by the user beforehand. For each heater voltage, v^, a non- 
linear equation, as described in Eq. (fl3 can be formed. By stacking N equations, a set of nonlinear 
equations is obtained as 



Y = 



vi 

V2 



PxS(X x ,vi) + m 

PxS{i x ^ 2 ) + m 

P x S{X Xl v N ) + r\ N 



(3) 



The unknown wavelength, X x , and power, P x , can be estimated by solving Eq. (0. Obtaining an 
analytical solution of Eq.© may be infeasible in practice. Alternatively, a numerical solution can be 
obtained by solving a nonlinear least squares (NLLS) problem as 



(l x ,P x ) 



arg min —\Y 



(4) 



with Y defined in Eq. (0 and Y defined as 



Y = 



h 



P.AX X7 V! ) 

P X S(k,V2) 

P x S(X x ,v N ) 



(5) 



Here S(X,Vk) is the approximated spectral sensitivity function that was obtained by calibration in 
Step 1 . The mean-square-error between the measurement vector Y and the approximation Y at certain 
guesses, X x and P x , is defined as the cost function / = jj\Y — Y\ 2 . Eq. (0]i aims to find an appropri- 
ate X x and P x such that the approximation Y is as close as possible to the measurement Y, i.e., J is 
minimized. In this way, the unknown wavelength, X x and power P x , of the incident light field are esti- 
mated, i.e., the wavelength readout is obtained. By accurately approximating the spectral sensitivity, 
S(X , Vk), and increasing the number of tunings (i.e., more equations are used), the accuracy of X x and 
P x can be improved. 

The special feature of our method is that the spectral sensitivity function, 5(A,v^), is not necessarily 
bijective (i.e., strictly increasing or decreasing) with respect to wavelength, X, in the operational range, 
as compared with the LUT method. Because of this, not only is the operational range of our wavemeter 
extended (non-monotonous regions are also allowed), but also the accuracy of the estimates is improved. 
Since monotonicity is not required, sensitivity curves with local peaks and troughs (and even multi-peaked 
and multi-troughed sensitivity functions) can be used for wavelength estimate. As more equations are added 
in Eq. (0, the accuracy of the estimate improves. 

3. Experiments and results 

For an experimental demonstration of our wavelength estimation method, the available MRR (with a radius 
of R=85.4 /im) is connected to an independently calibrated, fiber-coupled tunable laser source (HP8 1689A, 
Agilent). The light transmitted through the MRR is coupled directly to a photodiode (FGA10, Thorlabs). 
The signal of the photodiode (i.e., y^) is amplified and recorded by a computer via a data acquisition 
card (PCIe-6251, National Instruments). The heater voltage vj is provided by the computer and a power 
amplifier. 

To calibrate the sensitivity function, S(X,v), and to test the wavelength estimate method, two data sets 
were collected, denoted as {X,v,y c } and {X, v,y t }, respectively. These two data sets were generated using 
a set of known wavelengths and heater voltages. The laser was tuned over the wavelength range from 
1549 nm to 1553 nm, in 400 equal steps of 0.01 nm. The power was constant at 2 mW. For each wavelength, 
N = 21 heater voltages were applied, fromO V to 10 V in steps of 0.5 V, and the corresponding photodiode 
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Fig. 4. Variance accounted for (VAF) of the neural networks with respect to the number of 
neurons (Q). The VAF in the vertical axis is averaged over VAF of all 21 neural networks. 
As more neurons are used, the accuracy of the neural network increases until Q = 14. 



signals were collected twice, one for calibration (i.e., y c ) and the other for test (i.e., y t ). Thus, in total, 
401 x 21 = 8421 data points were recorded for calibration and the other 8421 data points for testing. A 
typical time required for collecting both data sets is about two hours, which is limited by the wavelength- 
tuning speed of the tunable laser. 

It can be seen in Fig. [3] that S(X , Vk) consists of a broad range of nonlinear variations vs. wavelength, 
resembling a set of mutually shifted Fabry-Perot resonator transmission functions. At longer wavelengths, 
the depth of the modulation decreases. This can be addressed to the overall heating of the optical waveguide 
chip, which causes the alignment between the optical filter and the waveguide to change over the course of 
the measurements. 

During calibration, each sensitivity function at a particular voltage as a function of wavelength is ap- 
proximated by a neural network, i.e., 21 neuronal network fits are performed, one for each heater voltage. 
To select an appropriate number of neurons, Q, in the network as presented by Eq. (O (2 layers, with Q 
neurons in the first layer and one neuron in the second layer), the variance accounted/or (VAF) was used 
as a criterion ll22l . The VAF is defined by 



A(y k ,yk) = i - 



\ ar(y k ~y k ) 
var(y k ) 



x 100%. 



(6) 



Here, A (%,yk) is the VAF between % and 3^, ranging from — °° to 100%. var(j^) is the variance ofyk- Fig.|4] 
shows the VAF of neural networks with respect to Q. As more neurons are used, the accuracy of the neural 
network increases until Q = 14. For Q > 14, the fitting accuracy is limited by the measurement noise in 
the data. Accordingly we decided to use 14 neurons in our network. The corresponding VAF (averaged 
over all 21 neural networks) is 99.63%, indicating an accurate approximation of the spectral sensitivity 
functions. This can be seen in Fig. [3j which shows an excellent agreement between the neural network 
fitting functions and the experimental data. 

After calibration, the wavelength and power estimate algorithm was tested with the data set {X,v,yt}, 
where X is considered unknown. For each "unknown" wavelength X x j, i e [1,401], the corresponding 
photodiode signals y t j.k,k = !,■■■ ,21 are scaled by a random number p X i € [0.1, 10]. This simulates an 
input that has both an unknown wavelength and an unknown power, which are to be estimated. The set 
of scaled photodiode signals are given as y t /# = pxfltikik = 1)"' ,21. Heating voltages, Vk, and signals, 
$t i k> were then fed into the wavelength estimation algorithm and the unknown wavelength A x .,-, and power- 
scaling factor p X i were estimated. 

Fig. [5] shows the wavelength estimation error X x — X x with respect to wavelength, X x , and Fig. [6] shows 
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Fig. 5. Wavelength estimation error X x — X x in the test set. 
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Fig. 6. Histogram of the wavelength estimation error. The maximal estimation error for the 
test set is 22 pm. 95% of the wavelength estimation errors are below 8 pm in the range 
between 1549 nm and 1553 nm. 



the estimation error in a histogram of occurrence. It can be seen that the maximum absolute wavelength 
estimation error for the test set is 22 pm and that 95% of the wavelength estimation errors are below 8 pm 
in the range from 1549 nm to 1553 nm. Fig. |7]shows the histogram of the relative power estimation error 
(px — Px)/Px with p x the estimate of p x . The error in power estimation is 1 .6% with respect to the known 
value; and in 95% cases, the power estimation error is less than 0.8% with respect to the known value. 

To analyze whether the proposed neural network and nonlinear least square (NN+NLLS) method pro- 
vides improved accuracy over a simple LUT method, we chose a LUT where the ratio between two spectral 
sensitivity curves strictly decreases. For the MRR used here, we found a spectral region that satisfied this 
constraint, e.g., between 1550.3 nm and 1550.8 nm for the heater voltages v^=0 V and 5 V. The wavelength 
was estimated using the LUT method and for the same data using the NN+NLLS method, as described 
above. The methods are compared in Fig. [8] The maximum wavelength estimation error in the named range 
is 42 pm for the LUT and 6 pm for our NN+NLLS approach. 95% of the errors are less than 27 pm for the 
LUT and less than 5 pm for NN+NLLS. This shows the NN+NLLS method, indeed, provides improved 
accuracy over that of LUT. The latter corresponds to a relative precision of about 3- 10 -6 . The improvement 
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Fig. 7. Histogram of the relative power estimation error (p x — Px)/px in the test data set. 
The relative power estimation error is at maximum 1.6%, and in 95% of the cases less than 
0.8%. 
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Fig. 8. Comparison between LUT and the proposed NN+NLLS method. 95% of the errors 
fall below 27 pm for LUT and 5 pm for the proposed method (about one fifth of that for 
LUT). 



is as big as a factor of five. 



4. Conclusions 

In conclusion, we have presented a novel method to measure the unknown power and wavelength of a laser 
with high precision and an extended measurement range. The method has been verified in a simple setup, 
based on an integrated optics micro ring resonator, which allows the wavemeter be integrated with other on- 
chip optical components. Neural networks were used to approximate the spectral sensitivity functions of the 
optical channel at different heating voltages. When injecting light with an unknown wavelength and power, 
nonlinear equations are formed and solved to provide an estimate of the wavelength and power. In our 
experimental verification, we demonstrated a spectral precision of about 8 pm (5-10~ 6 ) in the wavelength 
range between 1549 nm and 1553 nm. A higher precision of 5 pm (3-10~ 6 ) was achieved in the range 
between 1550.3 nm and 1550.8 nm, which is a factor of five compared to a simple lookup of data. 



The advantage of our approach is that it does not attempt to obtain a high precision and an increased 
wavelength range with a sophisticated and highly precise design of an appropriate spectral response func- 
tion. Instead, our approach is able to handle all deviations that may enter via fabrication processes or fiber 
input and output coupling, as long as the optical transmission function is reproducible. Essentially, the ap- 
proach replaces the need for a complex and precise optical design with a smart readout ("trading hardware 
for software"). This can help to reduce the cost of fabrication and expand the range of applications. 



